Cooperative domain type interlayer sp 3 -bond formation in graphite 
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Using the classical molecular dynamics and the semiempirical Brenner's potential, we theoretically 
study the interlayer a bond formation, as cooperative and nonlinear phenomena induced by visible 
light excitations of a graphite crystal. We have found several cases, wherein the excitations of certain 
lattice sites result in new interlayer bonds even at non-excited sites. We have also found that, a new 
interlayer bond is easier to be formed around a bond, if it is already existing. As many more sites 
are going to be excited, the number of interlayer bonds increases nonlinear ly with the number of 
excited sites. This nonlinearity shows 1.7 power of the total number of excited sites, corresponding 
to about three- or four-photon process. 

PACS numbers: 64.70.Nd,61.48.Gh,31.15.xv,05.40.-a 
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I. INTRODUCTION 

Recently, the phase transition from the sp 2 -bonded 
graphite to the sp 3 -bonded diamond induced by visible 
light irradiation has been experimentally proposed, and 
has attracted considerable attention. 1 As well known, 
the diamond is synthesized macroscopically from the 
graphite under high temperature and pressure (3000 °C, 
15 GPa),-^ or the irradiation of strong x ray 4^, due to 
the high energy barrier between them. It is estimated to 
be about 0.3-0.4 eV/atom by the local-density approxi- 
mation (LDA). 6,7 However, Kanasaki et at. found that, 
after the irradiation of visible photons onto graphite, a 
new nanoscale domain with sp 3 structure, including more 
than 1000 carbons, appears locally in the graphite surface 
and stable for more than 10 days at room temperature^ 
The nanoscale domain is an intermediate state between 
the graphite and the diamond. Therefore, it is called 
"diaphite i^" Further pulse excitations will generate a 
lot of such domains, and they combine with each other, 
proliferate, and are expected to result in the macroscopic 
diamond phase finally. In general, this kind of process is 
called the photoinduced structural phase transition. 10 

The essential aspects of this experiment are as follows: 
(1) the exciting laser, with an energy of 1.57 eV, should 
be polarized perpendicular to the graphite layer, while 
the one polarized parallel to the layer gives no effect. It 
means that only the interlayer charge-transfer excitation 
can trigger this process. (2) The exciting light should 
be a femtosecond pulse, while the picosecond one gives 
almost no contribution. It means that only a transient 
generation of an excited electronic wave packet in the 
semimetallic continuum can efficiently trigger this pro- 
cess. (3) The process is quite nonlinear but less than the 
ten-photon process. 

From the above experiments, we can think of the fol- 
lowing scenario for the present process £ when the fem- 
tosecond pulse is irradiated onto the graphite perpendic- 
ularly, an electron-hole pair spanning two layers is gen- 
erated. This electron-hole pair mainly dissipates into the 
semimetallic continuum of the graphite as plus and minus 



carriers due to the good conductivity of graphite. How- 
ever, by a small but finite probability, this electron- hole 
pair is expected to be bounded with each other through 
the interlayer Coulomb attraction. This exciton-like state 
self- localizes at a certain point of the layer by contracting 
the interlayer distance only around it. As the local con- 
traction of the interlayer distance proceeds, an interlayer 
a bond is formed. Through further pulse excitations, a 
lot of interlayer a bonds are formed stepwise, and then 
the diaphite structure is expected to appear macroscop- 
ically. 

In order to clarify the mechanism of this photoin- 
duced phase transition from the graphite to the diaphite, 
we have already calculated the adiabatic path from the 
graphite to the small diaphite domain, by means of the 
LDA (Ref. 11) and also by the semiempirical Brenner's 
theory ^ 12 i 13 These results have already well clarified the 
adiabatic property of the phase transition. We have also 
studied the early stage dynamics of the interlayer a bond 
formation along with the above scenario j 13 i 14 In the pre- 
vious paper ^ by solving the quantum time development 
of an electron-hole-phonon system, we have described 
the self-localization of an photoexcited electron-hole pair 
spanning two layers, full quantum mechanically. Since 
the subsequent interlayer a bond formation after the self- 
localization can be taken into account in a classical pic- 
ture due to the translational symmetry breaking of the 
excited stated we can calculate the bond formation us- 
ing classical molecular dynamics (MD) with the Bren- 
ner's potential. Consequently, we have found that the 
self-localization occurs by the probability of about 2% 
when the electron-hole pair is excited as a transient state 
with the energy of 3.3 ± 1.8 eV, while the subsequent 
bond formation is achieved when the excitation energy is 
more than 4.5 eV, corresponding to about three visible 
photons with the energy of 1.57 eV. 

So far, we have already clarified the dynamics how an 
interlayer bond is formed. The next stage is to investigate 
how a lot of interlayer bonds are cooperatively formed by 
further pulse excitations, that is, the effect of multisite 
excitation. In the case where several sites in graphite 



layers are excited, various cooperative phenomena will 
be expected to influence the interlayer bond formation. 
To confirm it, we calculate the dynamics of the multisite 
excitation. 

In the starting point of our calculation, we assume that 
the self-localization has already completed at more than 
one site in the graphite layers by the irradiation of many 
visible photons. Hence, in this paper, all the calculations 
of the bond formation are performed by using the MD, 
just same as in the previous paper. 14 Similarly, we adopt 
the Brenner's potential I in its original version, 15 which 
can well describe various carbon cluster systems. 



II. MODEL AND METHOD 
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FIG. 1: Schematic representation of AB stacking of two 
graphite layers. The solid line indicates the first layer and 
the dashed line is the second one. The carbon pair in which 
two carbons locate at the same site in xy plane is referenced as 
an a pair, and the pair which has a shift between two carbons 
in x direction is referenced as a f3 pair. 



We consider only two graphite layers with AB stacking 
and assume that the interlayer distance contracts locally 
at several sites of the crystal due to the self-localization 
of electron-hole pairs. Then, the two opposite carbons 
at each excited site are intruded inside of the two layers 
and have a velocity toward the inside. As an initial con- 
dition, therefore, we give the carbon pair initial intrusion 
(potential energy) and an initial kinetic energy. For sim- 
plicity, we also assume that the initial displacements and 
the initial velocities of the two carbons of the pair have 
the same magnitude in opposite directions. The initial 
kinetic energy Kq is given by Kq = 2^m c ^o, where m c is 
the mass of a carbon atom and vo is the initial velocity. 
The initial potential energy is obtained from the relation 
between the intrusion and the potential energy, calcu- 
lated by using the Brenner's theory ^ In our calculation, 
the total energy is always conserved, since we do not take 
into account any energy dissipation effects. The system 
consists of two graphite layers with the initial interlayer 
distance 3.35 A, and each layer includes 6240 carbons 
in about 128 A x 128 A, wherein a periodic boundary 
condition is imposed. 

First, we perform the MD calculation in the case where 
two or more carbon pairs are excited simultaneously in 
order to confirm the cooperative phenomena between ex- 
citations. In this calculation, we give the same potential 
energy and kinetic energy to all the excited carbon pairs, 
as an initial state. We calculate the dynamics for various 
spatial arrangements, various numbers of excited sites, 
and various excitation energies. 

Second, we investigate the easiness of interlayer bond- 
ing when an old interlayer bond has already existed 
nearby. As an initial state, we give excitation energy 
to a carbon pair near a certain interlayer bond, which is 
photogenerated beforehand. In this sense, the situation 
is a kind of multisite excitation, although the time inter- 
val between these two excitations is infinite. We perform 
the MD calculation for various values of excitation en- 
ergy and various positions of the excited carbon pairs. 
Thereby, we can estimate the lowest energy necessary to 
form an interlayer bond as a function of the distance be- 
tween the two excited sites. 



Finally, we consider the many and random multisite 
excitations whose calculation procedure is as follows: (1) 
we choose N ex excited sites randomly according to ex- 
citation density. (2) we determine the excitation energy 
E ex of each site randomly according to a Gaussian distri- 
bution G(E m ,AE) = exp[-{(E ex - E m )/AE} 2 }, where 
E m and AE are the mean energy and the width, and are 
set to be 3.3 eV and 1.8 eV, respectively. These values 
correspond to the energy of the transient state of the 
photoexcited electron-hole pair spanning two graphite 
layers, and have been obtained in our previous paper J^ 
We assume that this randomness of the excitation energy 
originates from the uncertainty of the transient excited 
state. (3) we determine the initial potential energy P and 
the kinetic energy K randomly under the condition that 
P + K = E ex . (4) we calculate the MD, starting from the 
initial state determined above, and count out the number 
of interlayer bonds formed after the time 1.0 ps. Within 
this time, the interlayer bond formation is expected to be 
almost complete. (5) we performed a similar calculation 
50 times and average them. 

In the procedure (1), due to the AB stacking, there 
are two types of carbon pairs, as shown in Fig. [TJ One 
is the case that the positions of two carbons at the first 
layer and at the second one are mutually same (a pair). 
Another is that the positions shift each other in x direc- 
tion (/? pair). In the latter case, the initial intrusion and 
velocity of the two excited carbons have not only the z 
component but also the x component, as indicated by 
arrows in Fig. [TJ Assuming only the x polarized light 
irradiation, all the j3 pairs have the shift only in the x 
direction. 



III. RESULTS AND DISCUSSIONS 

A. cooperative phenomena 

We show the results for the cases where two or more 
carbon pairs are excited simultaneously. From the results 
of the calculations for various initial states, we have found 
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FIG. 2: The example of cooperative phenomena in the case 
of 3 excited alpha pairs with the energy of 4.22 eV/pair. (a) 
Schematic representation of two graphite layers. The a pairs 
indicated by red closed circle are excited initially, having the 
kinetic energy of 4.08 eV/pair and the potential energy of 
0.14 eV/pair (the intrusion of 0.12 A), (b) The bond length 
between two opposite carbons as a function of time. The 
black solid line indicates the excite pair and the red dashed 
line is the central pair marked by a purple open circle in (a), 
(c) The snapshots of the MD calculation at time t — and 
0.5 ps. 



various interesting cases. Here, we discuss three examples 
for such cases shown in Figures [2j [3] and HI The first 
example is the case where three a pairs indicated by red 
closed circles in Fig. [2fa) are excited. Each pair has 
the initial kinetic energy of 4.08 eV and the potential 
energy of 0.14 eV corresponding to the intrusion of 0.12 
A. Hence, the total energy (4.22 eV/pair) is not enough 
to form an interlayer bond individually. We should note 
that the lowest energy of the bond formation is about 
4.5 eV, as mentioned above. That is, no interlayer bond 
can be formed in this energy region without cooperative 
effects of multisite excitation. 

Figure EJb) shows the time development of the bond 
length between two opposite carbons, where the black 
solid line indicates the originally excited pair and the 
red dashed line is the central pair marked by a purple 
open circle in Fig. EJa). As seen from this figure, the 
interlayer distance of the originally excited pair shrinks 
less than 2 A within 0.05 ps, resulting in an interlayer 
bond stable after 0.3 ps. With a little delay, even the 
interlayer distance of the central pair also shrinks, being 
dragged by the shrinkage of the bond length of the excited 
pairs, and results in a new interlayer bond. 

We show the snapshots of the MD calculation at t — 
ps and 0.5 ps in Fig. Efc). In this example, due to the 
cooperation of the three closely excited pairs, three inter- 
layer bonds can be formed at the originally excited sites, 
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FIG. 3: The example of cooperative phenomena in the case of 
3 excited alpha pairs with the energy of 2.51 eV/pair. The ki- 
netic and potential energies are 1.04 eV/pair and 1.47 eV/pair 
(the intrusion of 0.39 A), (a) The bond length as a function of 
time, plotted in the same way as Fig[2jb). (b) The snapshots 
of the MD calculation at time t — 0.5 ps. 



even though each excited pair does not have the suffi- 
cient energy to form the bond individually. Moreover, 
a new interlayer bond can be also formed at the central 
site, though it is not originally excited. This is because 
the extra energies from the originally excited pairs have 
propagated and concentrated at the center as a vibra- 
tional energy. This result means the nonlinearity that 
the number of interlayer bonds proliferates from to 4, 
because of the cooperativity. 

The second example is shown in Fig. [3j where three a 
pairs are excited similar to the previous example, while 
the initial kinetic and potential energies are 1.04 eV/pair 
and 1.47 eV/pair corresponding to the intrusion of 0.39 
A. Hence, the total energy is 2.51 eV/pair, being much 
lower than in the previous case. We also show the time 
development of the bond length and the snapshot at t = 
0.5 ps in Figures [3fa) and (b), respectively. As seen from 
these figures, the bond length of the excited pair shrinks 
a little but does not reach up to 2 A. Of course, this is 
because the excitation energy is too small to form an in- 
terlayer bond at the excited sites. However, the energies 
of the excited pairs concentrate on the central site and 
are superposed. Therefore, the interlayer bond is formed 
there at about 0.1 ps, as the bond lenth of the central 
pair shows. This example also shows the nonlinearity 
that the number of interlayer bonds proliferates from 
to 1. 

Similarly, we show the third example in Fig. [4j where 
six a pairs are excited. Compared with the previous ex- 
amples, the number of initially excited sites is larger, but 
each excited site is more distant from others, instead. 
The excitation energy is 4.74 eV/pair, in which the ki- 
netic and the potential energies are 3.20 eV/pair and 1.54 
eV/pair (the intrusion of 0.4 A). In Fig. HJ^b), the one- 
dotted-dashed line (blue) and the two-dotted-dashed line 
(green) denote the bond lengths of the intermediate a 
pairs numbered by 1 and 2 in Fig. HJa), respectively. As 
seen from Fig. HJ^b), the bond length of the excited pair 
shrinks rapidly and an interlayer bond is formed there, 
since the original excitation energy is enough to form this 
bond. Furthermore, the extra energies of these excited 
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FIG. 5: (a) The blue closed circle denotes the already existing 
interlayer bond and the excited sitesa pairs are numbered 
from 1 to 7. (b) the lowest limit of the excitation energy as 
a function of the distance between the already existing bond 
and the excited site. The horizontal dashed line indicates the 
lowest limit energy to form an interlayer bond individually. 



FIG. 4: The example of cooperative phenomena in the case of 
6 excited alpha pairs with the energy of 4.74 eV/pair. The ki- 
netic and potential energies are 3.20 eV/pair and 1.54 eV/pair 
(the intrusion of 0.4 A), (a) Schematic representation of two 
graphite layers, showing the initially excited a pairs indicated 
by red closed circle, (b) The bond length as a function of time, 
plotted in the same way as in the previous examples, while 
the one-dotted-dashed line (blue) and two-dotted-dashed line 
(green) denote the bond lengths of the intermediate a pairs 
numbered by 1 and 2 in (a), (c) The snapshots of the MD 
calculation at time t = and 0.5 ps. 



pairs propagate as far as to the central site gradually, 
through the intermediate sites as a vibrational energy. 
Then, the bond length of the central pair expands once 
but afterwards shrinks rapidly, resulting in the interlayer 
bond at about 0.2 ps. The nonlinear ity in this example 
indicates that the number of interlayer bonds proliferates 
from 6 to 7. 

Next, we discuss a more different cooperative phenom- 
ena from the previous ones. In the initial state of this 
case, an interlayer bond already exists at the site indi- 
cated by a blue closed circle in Fig. [5fa), and one of a 
pairs numbered from 1 to 7 by red is further excited. We 
performed the MD calculation starting from such an ini- 
tial state for various excitation energies and obtain the 
lowest excitation energy necessary to form an interlayer 
bond. Figure |5Jb) shows the lowest excitation energy 
as a function of the distance between the initially ex- 
isting bond and the newly excited site. The numbering 
by red corresponds to that in Figure [5fa). The hori- 



zontal dashed line at about 4.5 eV in this figure means 
the individual lowest energy of interlayer bonding. As 
seen from this figure, the shorter the distance from the 
existing interlayer bond is, the smaller the necessary ex- 
citation energy becomes. That is, an interlayer bond is 
easier to be formed around the existing bond. This is 
because carbons around the interlayer bond are intruded 
inside of the layers. When the excited pair locates about 
10 A away from the existing bond, the lowest energy is 
almost same as the case of the individual bond formation. 
These results also show the nonlinearity of the multisite 
excitation. 



B. random multisite excitation 

We have thus seen the several examples for coopera- 
tive phenomena. Next, we discuss the more realistic and 
random multisite excitation whose methodology has been 
already explained in detail in Sec. HU We calculate the 
random multisite excitation for several excitation densi- 
ties and show the results. In Fig. [6j as a example, we 
show the snapshots of the MD calculation for the excita- 
tion density 5 %, that is, the number of excited pairs is 
312, including a and /3 pairs randomly. These snapshots 
are taken from the view over the graphite layers. Gazing 
at them, a lot of six-membered rings in the layers can 
be recognized. In the left figures, the distance between 
two graphite layers is represented by the color palette 
whose relation with the distance is plotted at the left 
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FIG. 6: The snapshots of the MD calculation for the excita- 
tion density 5 % at t = 0, 0.33, 0.67 and 1.0 ps. The snapshots 
at the left side represent the distance between two graphite 
layers by the color palette at the left bottom, wherein the red 
means that an interlayer bond is formed. Those at the right 
side represent the shear of the upper layer to the lower layer 
and its direction by the circular color palette at the right bot- 
tom, wherein the white at the center of the circle means no 
shear and the shear of 0.71 A in the right direction is colored 
in red. 



bottom. Here, the red means that the distance becomes 
less than 2 A, that is, the interlayer bond is formed. In 
the right figures, the shear displacements of the upper 
layer relative to the lower one, together with its direc- 
tion, are represented by the circular color palette plotted 
at the right bottom. The white at the center of the circle 
means no shear. For example, if the carbon of the upper 
layer shifts 0.71 A relatively to the carbon of the lower 
one in x direction, these carbons are colored in red. 

Let us see the snapshots at the left side of Fig. [6] in its 
time order. At t — ps, the excited pairs are colored in 
yellow or orange since they are given some intrusion to 
the inside of the layers initially. At 0.33 ps, due to the 
motion of the excited pairs to the inside, the shrinkage 
of the interlayer distance occurs in almost all the area 
and the interlayer bonds appear at some of the excited 
sites. Afterwards, the part of the layers which does not 
contain the interlayer bonds separates away toward the 
original interlayer distance, as seen from the snapshot at 
0.67 ps. Finally, the system becomes almost stable with 
several interlayer bonds at 1.0 ps, wherein the shrinkage 
of the interlayer distance remains around the interlayer 
bonds. Of course, each carbon still keeps moving as a 
small vibration, since there is no energy dissipation in 
the system. In this 5 % excitation case, as seen from 
the snapshot at 1.0 ps, the interlayer bonds, which fi- 
nally remain, are rather rare. Therefore, the cooperative 
phenomena do not occur so intensively. 

Next, let us see the snapshots of the shear at the right 
side of Fig. [6] Since there is no shear at t — ps, the 
layers are white except the excited /3 pairs. At a glance, 
they might not be recognized to be white. However, by 
gazing at the figure in detail, it can be seen that the 
carbons and the bonds are colored in white. When a 
f3 pair is excited, the upper carbon shifts toward right 
relatively to the lower carbon, according to the initially 
given potential energy. Therefore, the excited /? pairs are 
colored in light red. As seen from the snapshots at 0.33 
and 0.67 ps, most of the carbons in the layers becomes 
red gradually, that is, the shear toward the right direction 
occurs in all the area of the layers due to the motion of the 
excited /3 pairs. It almost remains even at 1.0 ps. If two 
carbons of a /3 pair are stably bonded, these carbons keep 
red color. However, any interlayer bonds are not actually 
formed in j3 pairs. As a result of the MD calculation after 
1.0 ps, it is confirmed that the shear disappears gradually, 
and all the area of the layers return to be white again at 
about 1.5 ps, though we do not show it here. 

Similarly, as a example, we show the snapshots of the 
MD calculation for the excitation density 10 %, corre- 
sponding to 624 pairs' excitation, in Fig. [71 As seen from 
the left figures, the shrinkage of the interlayer distance 
occurs in almost all the area and the interlayer bonds 
appear at some of the excited sites, as well as in the case 
of 5 % excitation. While this shrinkage appears faster 
and the number of the interlayer bonds is more abundant 
than in the previous case, because of the larger number 
of excitations. The shrinkage finally remains only around 
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FIG. 7: The snapshots of the MD calculation for the excita- 
tion density 10 %, represented in the same way as Fig. \6\ 



the inter layer bonds at 1.0 ps. Compared with the case 
of 5 % excitation, the interlayer bonds cluster or close 
together. The cooperative phenomena occur much more 
intensively in this case. 

As seen from the snapshots at the right side of Fig. [7| 
the shear in the right direction also occurs in all the area 
once, but afterwards it disappears gradually. The occur- 
rence of the shear and its disappearance are also faster 
than in the 5 % excitation case. Moreover, the layers 
shift in the opposite direction only a little, as shown by a 
lot of carbons colored in light blue at 1.0 ps. This means 
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FIG. 8: The number of interlayer a bonds as a function of 
the number of excited sites N ex - The black dotted dashed 
line indicates the case where the interlayer bond increases 
proportionally to N ex , whose slope is 0.048. The blue dashed 
line is a fitted curve proportional to (Nex) 1 ' 7 • 



that no interlayer bond is formed in fj pairs. 

In Fig. [HI we show the number of final interlayer bonds, 
counted at 1.0 ps and averaged over 50 calculations, as a 
function of the number of excited sites N ex . It is plotted 
by the red solid line. The dotted dashed line (black) indi- 
cates the case where the interlayer bonds would increase 
linearly proportional to N ex . This slope is obtained as 
the probability that an interlayer bond is formed when 
one carbon pair alone is excited randomly. As seen from 
this figure, the number of interlayer bonds for more than 
5 % excitation increases nonlinearly with respect to the 
number of excited sites, and can be fitted by the blue 
dashed curve ex (Nex) 1 ' 7 . This means that effectively 1.7 
excitations contribute to the interlayer bond formation. 
While there is no cooperative phenomenon at all for less 
than 3 % excitation. In our calculation, the mean energy 
of the excitation is 3.3 eV, corresponding to about 2.1 
photons with the energy of 1.57 eV. That is, 1.7 excita- 
tions is equivalent to 1.7x2.1 = 3.57 photons on average. 
Therefore, we conclude that about three or four photons 
are effectively involved in this process. 

Incidentally, we discuss the shear between graphite lay- 
ers. The existence of the shear is very important for the 
diaphite domain formation, as discussed in the previous 



works, 



,11-13 



since interlayer bonds in /3 pairs are nec- 



essary to build up the diaphite structure. In our present 
calculation, the shear once appears but disappears as 
time passes. We think it is because the excitation in this 
case is a pulse, or only at the very beginning. As long 
as the light is irradiated for a certain period of time, 
the shear can be induced by the fast part of this pe- 
riod, and also the subsequent excitation can accelerate 
this shear by the later part of the period before it disap- 
pears. Therefore, it is probable that interlayer bonds are 
successfully formed in f3 pairs, so that the shear remains 



7 



stably, resulting in the diaphite domain. The calculation 
which takes into account the time width of the excitation 
is our future work. 



IV. SUMMARY 



the random multisite excitation. As a result, the num- 
ber of formed interlayer bonds increases nonlinearly with 
respect to the number of excited sites. The nonlinear- 
ity shows 1.7 power of the number of excited sites. This 
result means that about three or four photons are effec- 
tively involved in the interlayer bond formation. 



We have thus investigated the nonlinearity in the in- 
terlayer a bond formation, using the classical molecular 
dynamics and the Brenner's potential. Due to the multi- 
site excitation, various cooperative phenomena between 
excited sites appear, so that the interlayer bonds, more 
than the number of excited sites, are formed, or making 
up for the lack of energy each other. In the situation 
where an interlayer bond already exists, it is easier to 
form a new interlayer bond nearby. Based on the under- 
standing of such cooperative phenomena, we have studied 
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